Radial-breathing-like phonon modes of double-walled carbon nanotubes 
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The radial-breathing-like phonon modes (RBLMs) of the double-walled carbon nanotubes are 
studied in a simple analytical model, in which the interaction force constants (FCs) can be obtained 
analytically from the continuous model. The RBLMs frequencies are obtained by solving the dy- 
namical matrix, and their relationship with the tube radii can be obtained analytically, offering a 
powerful experimental tool for determining precisely the radii of the multi-walled carbon nanotubes. 
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Right now, the double-walled carbon nanotubes 
(DWCNTs) become an important research subject af- 
ter they are successfully synthesized by catalytic chem- 
ical vapor depositioni'- '■^ and the thermal conversion 
of Geo encapsulated in single-walled carbon nanotubes 
(SWGNTs)! Many researches^iii^i^ii^iaiiaiiiiHas have 

shown that there exists a kind of radial motion in the 
DWGNTs similar to the radial breathing modes (RBMs) 
in the SWGNTs, called as the radial breathing-like modes 
(RBLMs). But their relationship with the tube radii 
is still an open question, which is the most important 
for understanding the experimental RBLMs. In some 
publicat ions^'^'^'^ 1 10,11,12 ^ ^j^g force constants are com- 
puted numerically, making the relationship between the 
RBLM frequencies and the tube radii unclear. Dobardzic 
et al. suggested a simple analytical modeli^, considering 
the tube walls as coupled oscillators, and the FG values 
of the coupled oscillators are determined by a fitting to 
the experimental B2g graphite phonon frequency of 127 
cm~^. With this model, they have successfully calculated 
the RBLM frequencies of the DWGNT. 

In this work, we have investigated further the RBLMs 
of the DWCNTs, using the simple analytical modeli^, 
and assuming the carbon atoms on the nanotube wall are 
distributed continuously, which has been used in a lot of 
system s Based on the continuum model in 
Ref. [1^, we have derived a set of formulae to describe the 
interaction between an infinite length tube and an atom 
in any place. With these formulae, the force constants 
between the inner and outer tubes can be obtained an- 
alytically, and furthermore, an accurate relationship be- 
tween the RBLM frequencies and the tube radii has been 
obtained. Our final results only involve a few basic pa- 
rameter values relying on the nanotubes structure and 
the parameters of Lennard- Jones (LJ) potential. 

Now, let us introduce the simple analytical model. For 
the RBLMs, the only degree of freedom is the motion 
along the radial direction of the tube, so we can simply 
consider the RBLMs as the motion of a harmonic oscilla- 
tor. A segment of SWNT can be described as a smooth 
cylinder with its radius and height represented by r and 
h. Then its Hamiltonian can be written as: 



(1) 



Here the superscript S means the SWNT. m = 
(27Trh) na-rric is the mass of the cylinder, Ua- is the mean 
surface density of carbon atoms on the tube wall, and rric 
is the mass of carbon atom, r is the time derivative of the 
radius, and ojrbm is the vibration frequency of the RBM 
mode. It is well known that the RBM frequency of the 
SWNT satisfies a simple relationship with the nanotube 
radius, 



i^RBM = D/ (2ro) , 



(2) 



where D keeps a constant for any kind of SWNTs and tq 
is the equilibrium radius of the SWNT. Thus we have 
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For the DWCNTs, the Hamiltonian can be further 
written asi^: 
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where H^^^ and Hf^ are the Hamiltonians of the outer 
and inner SWNTs, respectively, with their radii as ri and 
r2 . And Vio is the interaction potential between the inner 
and outer tubes. 

In this work, we assume the interaction between inner- 
and outer-tube atoms can be described by the Lennard- 
Jones(LJ) potential. 
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with parameters, e = 2.964 meV and a — 3.407 A^^-. 
Based on our previous work^^, when the tube height is 
large enough, the interaction potential can be writtend 
as 
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Vio = * Snrieua * Vkernel 

+Nout*^TTr2en^ *Vkernel, (6) 



A'ira = 2'iTr2n„h and iVout = 2'iTrin„h are the number of 
atoms of the inner and outer tube segment, respectively. 
k is an integer. 



where 



Vkernel — 
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For the DWCNTs, when inner tube radius r2 is larger 
than (ri - r2) « 3.4 A, the < ^ 



so we can use our 



(8) series expression^^ for the Eq. ([8]), and get 
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Now we can obtain the dynamical matrix analytically, 



Here, 
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where so = ^, Si = ^, and S2 = ^. do = ri - r2 is 
the distance between the inner and outer tubes. Above 
formulae are suitable for any tube-tube distance. 

Assuming dp — 3.4 A« cr, which is equal to the inter- 
layer spacing of graphite, we can further get 



21 



K2 



= Ki2 = 



y/moutm.n dridr2 



Co 



1 9 1 

-16 H .s? H s? 

20 ^ 20 ^ 



1 d^V,o 



Co 



I6 + 8S2 



41 
20 



41 . 
40 



(12b) 



(12c) 



i^ii 



1 d^Vo 
T'out dr^ 



Co 



41 o 41 . 

16 - 8si s? s? 

20 ^ 40 ^ 



(12a) 



Here Co = 



Choosing = irT^-rr—, with 



ac-c the bond length of the SWNT, being about 1.42 
A, e = 2.964 meV, and rric = 12.011 amu, we can obtain 
Co w 241 cm^^. In this work, the prefactor D in Eq. [2] 
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is always chosen as 225 cm~-'^nm. By diagonalizing the 
dynamical matrix of Eq. 10, the RBLM frequencies can 
be obtained. 

For example, the RBLM frequencies of (5, 5)@(10, 10) 
DWCNT are calculated using above Co and D. Then one 
can easily obtain ojlf = 173.9 cm^^ and ujhf = 340.8 
cm~^. And the corresponding RBM frequencies of (5, 5) 
and (10, 10) SWNTs obtained by Eq. [2]are ^(5,5) = 331.9 
cm~^ and W(io.io) = 165.9 cm~^. Thus the frequency 
shifts of the RBLMs in the DWCNT are Auj^'^ =ujlf- 
^(1040) = 8.0 cm"i and A[j|^'^ = luhf - = 8.9 

cm"^, respectively. 

To make a comparison, a first-principle calculation has 
been performed on the (5, 5)@(10, 10) DWNT to check 
the validity of our method. A supercell geometry was 
adopted so that the DWCNT is aligned in a hexagonal ar- 
ray with nearby DWCNT center distance of 25 A, which 
is found to be larger enough to prevent the tube-tube in- 
teractions. The K-points sampling in the reciprocal space 
is a uniform grid 1 x 1 x 24 along the nanotube axis {Z di- 
rection) in our calculations. After structure relaxation on 
both of the lattice constant along the tube axis and the 
atomic positions, the optimal structure is obtained when 
the residue forces acting on all the atoms were less than 
0.01 eV/A. The final dynamical matrix is constructed 
using cumulant force-constant (CFC) method^SiiSl. Be- 
cause of the advantage of the CFC method in dealing 
with the low frequency part of the phonon dispersion, 
only a small supercell (1x1x2) is sufficient to calcu- 
late accurately the phonon modes on the F point and 
the X point. Our ah initio calculations were performed 
using highly accurate projected augmented wave (PAW) 
method^^, implemented in the Vienna ah initio simula- 
tion package (VASP) package^'^. They are based on the 
density-functional theory in the local-density approxima- 
tion (LDA). 

The obtained phonon dispersion curves and the den- 
sity of states for the (5, 5)@(10, 10) DWCNT are shown 
in Fig. 1, in which the low-frequency dispersions are 
very well reproduced. Its RBLM frequencies can be de- 
termined by analysis of the normal modes, which are 
ujfp = 182.2 cm-i and w^V = 354.6 craT^ . The 
RBM frequencies of (5, 5) and (10, 10) SWNTs are 
also calculated by the same ah initio method under 
the same calculation conditions as the DWCNT, which 
are w^g^'^^ = 343.8 cmT^ and wj'i^o.io) = l''3.6 cmT^ . 
Then, the corresponding frequency shifts of them in 
the DWCNT are Awf^ = ujfp - uj<^^q^^q) = 8.6 cm"! 

and Auj]fp — uj^fp — w^g^g^ = 10.8 cm^^, respectively. 
These results are consistent with our above obtained re- 
sults. The residual difference can be further minimized 
by adopting the first-principle tube radii of -R(5.5) = 3.411 
A and -R(io,io) = 6.766 A, and the first-principle RBM 
frequencies of the SWNTs in the dynamical matrix of 
Eq. 10. By this way, the new frequency shifts obtained 
by our analytical method will be Awf;^^ = 9.0 cm~^ and 
Alu^'p = 10.3 cm~^, agreeing very well with above first- 



FIG. 1: (Color online) Calculated phonon dispersion and the 
density of states for the (5, 5)@(10, 10) are shown in the left 
and right panels, respectively. 
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principle values. 

Our approach can be further applied to the multi- 
walled carbon nanotubes (MWCNTs). Here we take a 
3-layer MWCNT as an example. Again, we can write 
the system Hamiltonian as 



ff^ = E^f + E^«.' (13) 

1=1 



where Hf is the Hamiltonian of the ith layer, and Vij 
is the interaction potential between the ith and the jth 
layer nanotubes. And the radius of the ith layer is repre- 
sented by ri with ri > r2 > r^. So, its dynamical matrix 
can be written as: 
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where J^Tq?^ means the (a, /?) element of the K matrix, 
which comes from the interaction between the ith and 
the jth layer tubes. If condition 0^3 = l~^^J 1 can 
be satisfied, the Eq. (fTTj) can also be used to get the 
force constants between the 1st and the 3rd layer tubes. 
So using Eq. (fTTjl . the force constants between any two 
different layer tubes can be obtained easily. Here, for 
the 3 layers MWNT, we use do = 6.8 A as the distance 
between the 1st and the 3rd layer tubes, then Eq. ITU) 
is fully determined. 

To examine the precision of our result, we com- 
pare them with other numerical ones. The RBLM fre- 
quencies as a function of the outer tube radius for 
{n,n) @ {n + 5,n + 5) @{n + 10,n + 10) MWCNTs are 
plotted in Fig. 2 (For comparison, please see figure 5 
of Ref. [T3 ). From Fig. 2, it is found that our results 
match Ref. [12 very well, and the residual difference may 
come from the choice of some parameters, e.g., D, do 
and n^. In fact, if we only slightly adjust both of D 
and Co, our RBLM frequencies can match all the other 
publications^iii^iSiiaiiiii^ perfectly 



FIG. 2: (Color online) Dependence of the RBLM freqeun- 
cies on the outermost tube radius for 3 layer MWCNTs, i.e., 
(n, n) @ {n + 5,n + 5) @ (n + 10,n + 10). The highest (HF), 
the middle (MF) and the lowest (LF) RBLM are plotted. 
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Furthermore, we show the dependence of the RBLM 
freqeuncies on the number of layers in Fig. 3, in which 
the radius of the innermost tube is fixed at 30A. This 
result matches the Fig. 6 of Ref. [l^ perfectly. Based 
upon Fig. 3, a new method to determine the number 
of layers for the MWCNTs can be proposed, i.e., if the 
position of the lowest RBLM for a MWNT is obtained, 
one can easily find out how many layers are present in 
the MWCNT by comparing the experimental data with 
Fig. 3. 



FIG. 3: (Color online) Dependence of the RBLM freqeuncies 
on the layer number. Here, the radius of the innermost tube 
is fixed at 30 A. 



150 



a 100 



a 

a 
> 

03 

a 

a 
Pi 



50 



Out-of-plane counter-phase 
B^^ mode of bulk graphite 



• * • • 

: • : : 
- ■ ■ • 



m 



• 



6 8 10 
Number of layers 



12 



14 



Now we would like to illustrate the application of our 
approach. In Ref. i5|, the Raman investigation of the 
DWCNTs prepared from thermal conversion of Cgo en- 
capsulated in the SWCNTs has been reported, in which 
the Eq. (HD, suitable only to the SWCNTs, was used 
to obtain the tube radii. Now we used our approach to 
determine the inner tube radius and the tube-tube dis- 
tance. Here we choose a pair of experimental data from 
Ref. H, ujlf = 150 cm~^, and cohf = 267 cm^^. By Eq. 
(fTTj) . the high and low RBLM frequencies for different r2 
can be obtained at the certain dg. Then comparing them 
with the experimental frequencies, we can get two differ- 
ent inner tube radii. In a real system, these two inner 
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tube radii value must be the same, and so we have to ad- 
just the distance do between the inner and outer tubes. 
Finally, a most reliable inner tube radius can be found. 
In Fig. 4, we finally set do to be 3.558 A, making the 
cross points of low and high frequencies with the experi- 
mental frequencies lie on one vertical line, from which a 
reliable inner tube radius of about 4.453 A is obtained, 
i.e., R^n = 4.453 A, and so Rout = 8.011 A. 

FIG. 4: (Color online) The Raman frequency vs the inner 
tube radius, which is obtained after adjusting the distance do 
between the inner and outer walls to about 3.558 A. 



In summary, we presented an analytically solvable ap- 
proach for the radial breathing-like phonon modes of the 
MWNTs, where the tube walls are treated as coupled 
oscillators, and the carbon atoms on nanotubes are dis- 
tributed continuously. Then the force constants between 
inner and outer tubes, and the relationship of the RBLM 
frequencies with the tube radii can be obtained analyt- 
ically. It is found that our obtained RBLM frequencies 
are well consistent with those in other publications and 
experiments. Finally, our result can be helpful to deter- 
mine the radii of the DWCNTs and MWCNTs. 

Acknowledgments 



280 



260 



240 



s 



220 



3 

a 
> 

S3 



200 



180 



160 



140 




Most reliable 
inner tube radius 

R. =4.453 (A) 
in * ' 




co=267 (cm"'') 



4.35 4.4 4.45 4.5 

imer Tibe Eadius (A) 



4.55 



The authors acknowledge support from the Natural 
Science Foundation of China under Grants No. 10474035 
and No. A040108. The authors also thanks support to 
this work from a Grant for State Key Program of China 
through No. 2004CB619004. Our first-principle calcula- 
tions have been done on the Sun Fire V20z computers. 



* Electronic address: |E-mail: wu gaxp@gmail.com| 
t Corresponding author E-mail: jdong@nju.edu.cn 

^ T. Sugai, H. Omote, S. Bandow, N. Tanaka, and H. Shi- 

nohara, J. Chem. Phys. 112, 6000 (2000). 
^ E. Flahaut, A. Peigney, Ch. Laurent, A. Rousset, J. Mater. 

Chem. 10, 249 (2000). 
3 H. Zhu, C. Xu, B. Wei, and D. Wu, Carbon 40, 2021 

(2002). 

* S. Bandow, M. Takizawa, H. Hirahara, M.Yudasaka, and 
S. lijima, Chem. Phys. Lett. 337, 48 (2001). 

^ S. Bandow, G. Chen, G. U. Sumanasekera, R. Gupta, M. 
Yudasaka, S. lijima, and P. C. Eklund, Phys. Rev. B 66, 
075416 (2002). 



^ E. Dobardzic, I. Milosevic, T. Vukovic, B. Nikohc, and M. 

Damnjanovic, Eur. Phys. J. B 34, 409 (2003). 
^ M. Damnjanovic, I. Milosevic, E. Dobardzic, T. Vukovic, 

and B. Nikolic, Phys. Rev. B 69, 153401 (2004). 
* M. Damnjanovic, E. Dobardzic, 1. Milosevic, T. Vukovic, 

and B. Nikolic, New J. Phys. 5, 148 (2003). 
^ M. Xia, S. Zhang, E. Zhang, S. Zhao, and X. Zuo, Phys. 

Rev. B 69, 233407 (2004). 
^° R. Pfeifler, Ch. Kramberger, F. Simon, H. Kuzmany, V.N. 

Popov, and H. Kataura, Eur. Phys. J. B 42, 345 (2004). 

A. Rahmani, J.-L. Sauvajol, J. Cambedouzou, and C. 

Benoit, Phys. Rev. B 71, 125402 (2005). 

V. N. Popov and Luc Henrard, Phys. Rev. B 65, 235415 



6 



(2002). 

E. Dobardzic, J. Maultzsch, I. Milosevic, C. Thomsen, M. 
Damnjanovic, Phys. Status Solidi B 237 R7 (2003). 
L.A. Grrifalco, J. Phys. Chem 96, 858 (1992); L.A. Giri- 
falco, Phys. Rev. B 52, 9910 (1995). 

C. Rey, J. Garcia-Rodeja, L.J. Gallego, J. A. Alonso, Phys. 
Rev. B 55, 7190 (1997). 

J. Song and R.L. Cappelletti, Phys. Rev. B 50, 14 678 
(1994). 

L. Henrard, E. Hernandez, P. Bernier, and A. Rubio, Phys. 
Rev. B 60, R8521 (1999). 

L.A. Girifalco, Miroslav Hodak, and Roland S. Lee, Phys. 



Rev. B 62, 13104(2000). 

G. Wu and J. Dong, Phys. Rev. B 71, 115410 (2005). 
K. Parlinski, Z-.Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 
78, 4063 (1997). 

21 L.H. Ye, B.G. Liu, D.S. Wang, and R. Han, Phys. Rev. B 
69, 235409 (2004). 

22 RE. Blochl, Phys. Rev. B 50, 17 953 (1994); G. Kresse 
and D. Joubert, ibtd. 59, 1758 (1999). 

G. Kresse and J. Hafner, Phys. Rev. B 48, 13115 (1993); 
G. Kresse and J. Furthmiiller, Comput. Mater. Sci. 6, 15 
(1996); Phys. Rev. B 54, 11 169 (1996). 



